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ABSTRACT 

We study heating of cool cores in galaxy clusters by cosmic-ray (CR) stream- 
ing using numerical simulations. In this model, CRs are injected by the central 
active galactic nucleus (AGN) and move outward with Alfven waves. The waves 
are excited by the streaming itself and become non-linear. If magnetic fields 
are large enough, CRs can prevail in and heat the entire core because of a large 
Alfven velocity. We find that the CR streaming can stably heat both high and 
low temperature clusters for a long time without the assistance of thermal con- 
duction, and it can prevent the development of massive cooling flows. If there is 
even minor contribution of thermal conduction, the heating can be more stabi- 
lized. We discuss the reason of the stability and indicate that the CR pressure 
is insensitive to the change of intracluster medium (ICM), and that the density 
dependence of the heating term is similar to that of the radiative cooling. 

Subject headings: cooling flows — cosmic rays — galaxies: clusters: general - 
waves — X-rays: galaxies: clusters 



1. Introduction 



Clusters of galaxies are filled with hot X-ray gas (intracluster medium; ICM). Although 
the cooling time of the ICM is larger than th e age of the Universe for the most part of a clus- 
ter, the core is the exception (jSarazinlll986l ). Since the cooling time at the core (~ 10 8 yr) 
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is much smaller than the age of the cluster (> 10 9 yr), a subst antial gas inflow, which 



was called a "cooling flow", was expected to develop (|Fabianlll994j . and references therein) 



from intermediate or low temperature gas (e.g. 


Ikebe et al. 


1997 




Makishima et al. 


Peterson et al. 


2001; 


Tamura et al. 


2001; 


Kaastra et al. 


2001; 


Matsushita et al. 


2002h 



This 

means that the the cooling rate is much smaller tha n that previously assumed. Chandra ob 



servations are also consistent with these results (e.g. 


McNamara et al. 


2000; 


Johnstone et al. 


2002 




Ettori et al. 


2002 




Blanton, Sarazin. & McNamara 


2003 


). These observations show 



that some heat sources balance with radiative cooling and prevent the development of cool- 
ing flows. 

Many mechanisms for heat sources of the ICM have b een proposed. Thermal conduction 



from the hot outer laye rs of clusters is one popular idea (ITakahara &; Takaharal 1 19791 . 11981 



Tucker fc Rosnerlll983l ). This model may work for c 
ture if the conductivity is ~ 30% of t he Spitzer one ( 



the conduct i vity requires fine -tuning (IBregman fc David 



usters with a middle or hig h tempera 



Zak amska fc Naravanl 120031). However 



1988 



Soker 


2003: 


Guo & Oh 


2008) 



Brighenti & Mathewsl 12003 



creases for low temperature clusters, thermal conduction may not be able to transfer enough 



heat 


sources (e.g. Churazov et al. 2001 


Quilis. Bower. & Balogl 


2001; Briiggen & Kaisei 




2002; 


Basson & Alexander 


2003). In fact, it has been observed that the AGNs c 


isturb 


the surrounding ICM 


(e.g. 


Fabian et al. 


2000l McNamara et al. 


2000: iBlanton et al. 


2001 




McNamara et al.boOl 


: Mazzotta et al.lhood: Fuiita et al. 2002: Johnstone et al.ll2002: ] 


KemDner. Sarazin. & 


2002|; 


Takizawa et al. 


2003; 


Fujita et al. \ 


2004bf). However, it is not clear how the energv een- 



erated by the AGN is transfered to the ambient ICM. For example, sound waves observed 
around the AGNs in some clusters were considered the sign i ficant channel for the AGN 
energy input into the ICM (IFabian et al.l 120031 ; iForman et al.ll2005l ). However, theoretical 
studies have indicated that the dissipation of the sound waves or weak sho cks is too fast to 
heat an entire cool core and thus the waves ca nnot sta bly heat the core (jFuiita fc Suzuki 
20051 : iMathews. Faltenbacher. fc Brighenti! I2OO6I : see also iFuiita et~al"1l2007bh . 



(e.g. 


Tucker & Rosner 


1983; 


Bohringer & Morfil] 


1988; I 


lephaeL 


19871: 


ReDhaeli k, Silk 


1995; 


Colafrancesco, Dar, & De Riiiula 


2004; 


Pfrommer et al. 


20071 J 


ubelgas et al. 


2008|). In par- 



ticular, CR streaming has been investigated as a way through which the CR energy i s 
transfered to the ICM (IBohringer fc Morfilllll988l ; lLoewenstein. Zweibel. fc Begelmanlll99ll ). 
In this mechanism, PdV work done by the CRs on an Alfven wave effectively becomes an en- 
ergy source of the ICM (see §12.21). It is to be noted that this mechanism may also be wo rking 
in supernova remnants (e.g. IVolk. Drury. fc McKenzidll984l ; iBerezhko fc Ellison! Il999l ). For 
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clusters of galaxies, the models of iBohringer fc Morfilll (119881 ) and lLoewenstein et al.l (ll99lf ) 
are time- i ndependent . Thu s, the models cannot treat the ICM that is not in a steady state. 
Recently, iGuo fc Ohl (120081 ) studied non-steady models. They combined CR streaming and 
thermal conduction as heat sources, and showed that the ICM turns out to be in a steady 
state, and that a cooling flow is suppressed. In their models, CR streaming and thermal 
conduction seem to equally contribute to the heating of the ICM (their Figure 6). Since 
they constructed models for a middle-temperature cluster (A 2199), it is not certain whether 
the models c an be applied to low temper ature clusters in which thermal conduction is not 



effective (see iBrighenti &: Mathews! 120031 ) . Moreover, the conductivity of the ICM may be 



small even for hot clusters. The existence of cold fronts in some clusters means that at 
least in some regions in those cluste rs, the conductivity has to be much smaller than the 
Spitzer value ( lEttori fc Fabianl 120001 ) . Sim ultaneous heating by both the central AGN and 
thermal conduction had been proposed by iRuszkowski fc Begelmanl (120021 ). In their model, 
mechanical heating through bubble motion was considered for the AGN heating. Since the 
advantage of CR streaming over the mechanical heating was not clear in the models of 
Guo fc Ohl (120081 ). it may be useful to study the case where a cool core is heated only by 
CR streaming. 

In this study, we revisit the heating of the ICM by CR streaming. We consider non- 
steady ICM and investigate whether the CR streaming alone can stably heat the ICM without 
the assistance of thermal conduction. We emphasize that we do not intend to search models 
in which the heating is comp letely balanced with rad iative cooling. Since the age of clusters 
is finite (say ~ 5 Gyr; e.g. iKitayama &: Sutol 119961 ). it is sufficient to find solutions that 



are stable for that time. We t reat both h i gh an d low temperature clusters. We follow the 
growth of Alfven waves, which IGuo fc Ohl (120081 ) did not consider. We simulate for various 
parameters in order to find the nature of the CR heating, although our goal is not to compare 
the results with specific clusters in detail, because supplemental heating mechanisms other 
than CR streaming (e.g. shock heating) are likely to be effective in real clusters. We consider 
protons as CRs. 

The paper is organized as follows. In § 2, we explain our models for CR streaming and 
galaxy clusters. In § 3, we present the results for various parameters, and show that the 
CR heating is fairly stable. In § 4, we discuss the reason of the stability. § 5 is devoted to 
conclusions. 
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Models 



are 



2.1. Basic Equations 

For simplicity, we assume that the cluster is spherically symmetric. The flow equations 

0, (1) 

(2) 
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(4) 



where p is the gas density u is the gas velocity P g is the gas pressure, P c is the CR pressure, 
Pb is the magnetic pressure, G is the gravitational constant, M(r) is the gravitational mass 
within the radius r, k(T) = f c KoT 5 ^ 2 is the coefficient for thermal conduction and T is the 
temperature, n e is the electron density, A is the cooling function, H st is the heating by CR 
streaming, H co \\ is the heating by Coulomb and hadronic collisions, u is the CR transport 
velocity, D(p) is the diffusion coefficient for CRs averaged over the CR spectrum, r loss is 
the energy loss by Coulomb and hadronic collisions, and S c is the source term of CRs. 
Energy densities of the gas and the CRs are respectively defined as e g = P g /(j g — 1) and 
e c = Pc/ (7c - 1), where 7 9 = 5/3 and 7 C = 4/3. 

For thermal conductivity, Ko = 5 x 1CT 7 in cgs units corresponds to the classical Spitzer 
val ue, and f c is the ratio to tha t. The cooling function is based on the detailed calculations 
bv lSutherland fc Dopital Jl993h . 



n z e A = [Cx{k B T) a + C 2 {k B Tf + C 3 ]n,n, 



(5) 



where rii is the ion number density and the units for ksT are keV. For an average metallicity 
Z = 0.3 Z Q the constants in equation fl5j) are a = —1.7, (3 = 0.5, C\ = 8.6 x 10 -3 , C 2 = 
5.8 x 10 -2 , and C3 = 6.4 x 10~ 2 , and we can approximate riin e = 0.704(p/m^) 2 . The 
units of A are 10~ 22 ergs cm 3 (jRuszkowski fc Begelmanll2002l ). For Coulomb and hadronic 



collisions, we use H, 



coll 



ri r n e e r . and Y 



loss 



( c n e e c , where rj c 



2.63 x 10" 16 cm 3 s" 1 



and 



( c = 7.51 x 10 16 cm 3 s 1 (IGuo fc Ohll2008l ). Since their contributions are minor, the details 



of the collisions do not affect the results. 



-5 - 



2.2. CR Streaming and Heating 



Alfven waves, which scatter CRs, are driven by CR streaming bec ause CRs can give 
their momentum to the waves through resonance (streaming instability; ISkillina Il975l ; iBell 
19781 ). Since the CRs as a whole move with Alfven waves, the CR transport velocity in 
equation (j3J) is given by u = u + va, where va = B / \/Attp is the Alfven velocity for a 



magnetic field B (e.g. iKang & Jonesll2006l : iGuo & Ohll2008l : ICaprioli. Blasi. & Amatoll2009h . 

The wave energy Ua = SB 2 /(8n), where 5B is the magnetic field fluctuation, is amplified by 
the PdV work done by the CRs on Alfven waves: 



dt 



3Pr 



dr 



(6) 



(ILucek fc Belli 120001 ) . We solve this equation by setting Ua = at t = 0. 



After the wave energy increases to Ua ~ Um, where Um = B 2 /(8n) is the energy of 
the background magnetic field, the waves are expected to be damped by non-linear effects. 
Although the exact value of the maximum of Ua has not been known, we set it at Um- 
After the wave energy increases to the maximum, the waves are expected to heat ICM (e.g. 
Ohira et al.ll2009l ; iGargate et al.l 120101 ). Thus, we give the heating term in equation ([3]) by 

Tv A (7) 



H. 



Ht 



dr 



(jVolk et al.lll984t IKang fc Jonesll2006l ). We expect that T ~ 1 after the wave energy increases 
to Ua ~ Um- Thus, we simply give Y = Ua/Um for Ua < Um and T = 1 after Ua reaches 
Um- Since Ua rapidly increases at the cluster core fse e § I3.2p. the resul ts are not much 
dependent on the definition of T for Ua < Um- Note that IGuo fc Ohl (120081 ) assumed that T 
is always one. 

As the source of CRs, we primarily consider the AGN at the center of the c luster. The 
CRs may be supplied from bubbles observed in the central region of clusters ( IGuo &: Oh 
20081 ). Moreover, they may also be supplied by strong outbursts of the AGN. In this case, 
CRs ar e accelerated at the forward shock of a cocoon and they are distributed in a broad 



region (IFuiita et al. 



clusters ( McNamara et al. 



2007al). In fact, this kind of ou tbursts have been observed in several 

2005 aj]bT). Furthermore, c l uster mergers gen 



2005: Nulsen et al. 



erate turbulence around the core (IFujita. Matsumoto. fc Wadal l2004al ; iFujita et al.l 12005 



Ascasibar fc Markevitchll2006l ). and the turbulence could accelerate and provide CRs in the 



2003 



central region of the clusters fe.g.lOhno. Takizawa. fc Shibatall 2002 



Brunetti et al1l2004i Icassano fc Brunettill2005l : 



EnBlin et al. 



Fuiita. Takizawa. &: Sarazin 



201ll ; lBrunetti fc Lazarian 



201ll ). Although the supply of the CRs may be intermittent, we study co ntinuous CR iniec- 



tion as the time-average. The source term we adopt is similar to that of IGuo &: Ohl (120081 ) 
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and it is simply given by 



Sr. 



^AGN 



47r rf(ri/ro) v — r\ \r 



- a 



-(r/r ) 2 



-(r/n) 2 



where Lagn is the energy injection rate from the AGN. This means that CRs are mostly 



injected at Tq < r < r±. We adopt r = 20 kpc following iGuo fc Ohl (120081 ) based on 



the observations of bubbles. We also a dopt 7j = 150 kpc, which is the size of the radio 
minihalo observed in the Perseus cluster (IGitti. Brunetti. fc Settill2002l ) . iFujita et al.l ( ]2007af ) 
indicated that the ICM might be heated within that radius. The injection rate is Lagn = 
eMc 2 , where e is the parame ter, M is the inflo w rate of the gas toward the AGN, and c is 
the speed of light. Following IGuo fc Ohl (120081 ) , we assume v ~ 3. 



The diffusion of a CR particle depends on its energy, because its resonance with an 
Alfven wave depends on the gyro radius. This means that the diffusion coefficient D depends 
on the energy spectrum of CRs. However, the energy spectrum of CRs ejected from the AGN 
is unknown. Therefore, we give it as a simple function of the ICM density: 



D(p)=D (p/p c 



(9) 



where Dq and po axe respectively the values of D and p at r = and t = 0, and d is the 
parameter. We included the dependence on p, because we expect that the dif fusion coefficien t 
is reduced as the ICM is compressed and the magnetic fields are increased ( ]Mathewsll2009l ). 
We assume that the magnetic fields are B = fi (p/po) d > where B Q is the parameter and we 
take d = 2/3. 



2.3. Gravitational Matter and Gas Profile 



We consider two types of clusters. One is a relatively large cluster and the other is 
a small cluster. For the profile of the acceleration of gravity or the profile of gravitational 
matter, we adopt models constructed on observations of the Perseus cluster (large) the Virgo 
cluster (small). 



For the Perseus cluster, iMathews et al.l (120061 ) gave an analytical profile of the accelera- 
tion o f gravity, GM(r)/r 2 , which is constructed based on the observations by lChurazov et al. 
( 200 4]) . We adopt the profile including the contribution of the central galaxy (§ 2 in 



Mathews et al.l 12006). For t h e Vir go cluster, we use the density and temperature profiles 
obtained by iGhizzardi et al.l (120041 ) . From their results, we can construct the mass profile 
assuming that the ICM is almost in pressure equilibrium. 
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We assume that the ICM is initially isothermal. This is because we do not know the 
initial distribution of CRs. If we assume that the ICM is isothermal, u — 0, and P c = 
at t — 0, the ICM density is relatively low at the cluster center, and M is small when t 
is small. Then, as the ICM cools at the cluster center, M and the activity of the central 
AGN gradually increase, and the CRs injected by the AGN are accumulated in the core. 
On the other hand, if we adopt the observed current density and temperature profiles, the 
density at the cluster center is high and the cooling time is small. If we start calculations 
with u = and P c = at t = 0, M abruptly increases and CRs are injected in a very short 
time, which causes numerical instability. In reality, it is likely that AGN activities precede 
cluster formation. Thus, it is natural to assume that some amount of CRs had already been 
injected when the cool core was established. 

At t — 0, we assume that the temperature of the large cluster is 7 keV and that of 
the small cluster is 2.4 keV, which are the values in the outer region of the Perseus and the 
Virgo cluster, respectively (IChurazov et al.l l2004j ; iGhizzardi et al.ll2004j ). The ICM profiles 
are built so that the ICM is in pressure equilibrium in the given gravitational fields. The 
normalization of the density is determined so that the density in the outer region of the 
clusters is identical to the observed ones. The initial velocity of the ICM is u = 0. There 
are no CRs at t — 0. 



3. Results 



3.1. Numerical Methods 



The hydrodynamic part of the e quations is solved by a second-o rder advection upstream 



splitting method (AUSM) bas ed on iLiou fc Steffenl (119931 . see also IWada fc Norman! 12001 



Fujita. Suzuki, fc Wadall2004d ) . We use 300 unequally spaced meshes in the radial coordinate 
to cover a region with a radius of 1 Mpc. The inner boundary is set at r min = 5 kpc. We 
adopt inflow/outflow boundary conditions at the inner and outer radii. Model parameters 
are presented in Table [TJ For gravitational potential, 'P' refers to the Perseus type cluster, 
and 'V refers to the Virgo type cluster. 

Before we investigate heating by CR streaming, we study a pure cooling flow model for 
comparison. If there is no heating source (e = and f c = 0; Model LCF0), a cooling flow 
develops and reaches almost a steady state at t > 4 Gyr. Figure [1] shows the evolution of 
M. At t = 12 Gyr, M increases to 760 M Q yr" 1 . 
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3.2. CR Heating 

In this subsection, we consider our fiducial model with CR heating (Model LCRO). We 
did not include thermal conduction or f c = 0. The initial magnetic field and the diffusion 
coefficient at the cluster center are B = 10 fiG and D = 1 x 10 26 cm 2 s -1 , respectively. We 
use the same value of D Q for other models including CR heat ing. The diffus i on coe fficient is 



much smaller than the values assumed by iGuo fc Ohl (120081 ) and iMathewd ( 120091 ) . This is 



because the Alfven waves become non-linear (SB ~ B) as we show bellow. In th is case, the 



diffusion coefficient is close to the one for the Bohm diffusion (equation [14] in iBelll 119781 ) 
and can be very small. Effectively, the diffusion coefficient we adopted is too small to affect 
the results. In other words, the results are not much different even if we assume D = 0. 

Figure |2] shows the profiles of ICM temperature and density for Model LCRO. For 
t > 4 Gyr, they do not much change. However, the temperature of the inner boundary 
r = '"mm slowly oscillates at 0.6 < /c#T < 1 keV. The oscillation is reflected in M (Figured]). 
Compared with the pure cooling flow model (Model LCF0), M is significantly reduced (~ 
120 M yr" 1 at t ~ 12 Gyr), which means that CR streaming can be an effective, stable 
heating source. Figure [3] shows the evolution of the ratio Ua/Um- The region where Ua/Um 
reaches one expands inward and outward in the cluster. This means that the Alfven waves 
become non-linear in a wide region of the cluster at the end of the calculation. It is to be 
noted that even if we assume Ua/Um = 1 throughout the calculation, the results do not 
much change. The growth time of Ua is much larger than that for a supernova remnant. 
The main reason is that the spatial scale of a cluster is much larger than that of the precursor 
of the shock of a supernova remnant. The difference affects the gradient of P c in equation 

Figure H] shows the ratios P c /Pg and Ps/Pg at t = 9 Gyr. Both the CR and magnetic 
pressures are smaller than the gas pressure, although they are relatively large in the central 
region and cannot be ignored (P c /P g ~ 0.3 and Pb/P 9 ~ 0.3). In Figure |5j we show relative 
importance of the two heating mechanisms (H st and H co \\) at t = 9 Gyr. CR streaming alone 
can almost balance with radiative cooling except for the very inner region of the cluster. 
This means that it can heat almost the entire cool core of r < 100 kpc. The contribution 
of Coulomb and hadronic collisions are minor. Figure [6] shows the evolution of the ratio of 
the heating by CR streaming to the radiative cooling. The ratio gradually reaches one. The 
bend at r ~ 400 kpc at t — 6 Gyr corresponds to the point where Ua reaches Um (Figure E]). 
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3.3. Parameter Search 

In this subsection, we change model parameters to see how the results are affected by 
them. Probably, most uncertain parameters in our models are those for the energy input 
from the AGN. In Models LCRel and LCRe2, we change the value of e (Table [1]). The 
evolution of M is presented in Figure [7J As is expected, we tend to have a smaller M for a 
larger e. For Model LCRe2, the ICM becomes unstable for t > 8 Gyr, although M is smaller 
than that in Model LCRO for a long duration of ~ 5 Gyr (3 < t < 8 Gyr). The temperature 
and density profiles are shown in Figure [8] In our source model, CRs are injected most 
intensively at r ~ r (equation [8]). If e is too large, radiative cooling cannot cancel the CR 
heating at r ~ r , which makes the temperature and density profiles irregular at r ~ r . 
In general, models in which CRs are injected in a too narrow region and/or the ICM is too 
strongly heated tend to be unstable. Since our simulations are one-dimensional, we cannot 
investigate what happens after the ICM distribution becomes irregular. Mult i- dimensional 
simulations would be interesting to study that. The instability for Model LCRe2 can be 
prevented by thermal conduction. Model LCRc is the same as Model LCRe2 but f c = 0.1. 
In this model, the ICM is stable even at t ~ 12 Gyr. Moreover, M is smaller than that in 
Model LCRO throughout the calculation and ~ 40 M Q yr" 1 at t ~ 12 Gyr (Figured]). Note 
that if we do not include heat sources except for the thermal conduction of f c = 0.1, the 
mass inflow rate is M ~ 4OOM yr _1 at t ~ 12 Gyr. This means that the thermal conduction 
of this level alone cannot effectively halt a cooling flow. 

Models LCRnl and LCRn2 are the cases where the value of v is changed (Tabled]). The 
evolution of M is presented in Figure [9] For a larger z/, CRs are injected more intensively at 
r ~ r (equation [8J). Thus, the ICM becomes unstable at t ~ 9.6 Gyr for Model LCRn2. In 
Model LCRb, we change the strength of the background magnetic fields and adopt B = 5/iG. 
The results are not much different from those for Model LCRO (Bo = 10 /iG; Figure [9]). 
However, the ICM becomes unstable at t ~ 10.5 Gyr. 

We also considered a less massive cluster (the Virgo type). The initial ICM temperature 
(~ 2.4 keV) is much smaller than that of the Perseus type cluster (~ 7 keV). Model SCF0 
corresponds to a pure cooling flow. In Figure [101 the mass flow rate amounts to M > 
80 Mq yr -1 at the end of the calculation. Model SCR0 includes CR heating. The mass 
inflow rate is significantly reduced by the heating and M ~ 13 M & yr" 1 at t ~ 12 Gyr 
(Figure [TO]) . Figure [UJ shows the temperature and density distributions for Model SCR0. 
They are stable until the end of the calculation even if there is no assistance of thermal 
conduction. 
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Discussion 



The above results show that heating by CR streaming can almost balance with radiative 
cooling, and the heating process is relatively stable, even if there is no thermal conduction. 
The stability can roughly be explained as follows. 

In our calculations, the ICM velocity is much smaller than the Alfven velocity. Figure IT21 
is an example (Model LCRO at t — 9 Gyr). This means that va is the main contributer of 
u = u + v a in equation (J4]). Owing to the large Alfven velocity, CRs can prevail in and 
heat the entire core, which is different from other conventional heating me chanisms such as 
sound waves or weak shocks ( jFujita fc Suzukil 120051 ; iMathews et al.ll2006l ). Moreover, the 



-1/2 



cx p 



1 / 6 if we assume d = 2/3. Thus, the 



Alfven velocity is given by v a = B / ^Anp oc p d 
Alfven velocity is not much dependent on the ICM. The diffusion term in equation (j3J) can 
effectively be ignored because of the small value (§ 13.21) . These indicate that the distribution 
of P c or e c is insensitive to the change of local ICM, which makes the CR heating stable. 
Although the absolute value of P c increases as S c (oc M) increases, the overall shape of the 
profile of P c does not much change (Figure [T3]) . Moreover, since P c reflects accumulation of 
CRs injected so far, it is insensitive to a temporal change of S c . 

The stability also resides in the heating function H st = VA\dP c /dr \ and its global balance 
with the radiative cooling function. Here, physical quantities are the typical ones for < r. 
Since va does not much evolve and the overall shape of the profile of P c does not much change, 
the heating function can be approximated by H st oc Mf(r), where f(r) is a function of r and 
is almost independent of t. The mass inflow rate is given by M = 47rr 2 p\u\, which does not 
depend on the radius at a given time because of the mass conservation. The flow time of the 
ICM, t flow oc r/\u\, is nearly proportional to the cooling time of the ICM, £ coo i oc P g /{n 2 e h), 
because the flow compensates the cooled gas. Thus, we obtain H st oc (p 3 A / P g )r 3 f (r) . If 
the ICM is adiabatic. 



P g oc p 7s , where 7 9 



5/3. However, radiative cooling is effective in 
the central region of a cluster, and thus P g oc p 7 g, where 1 < iL < 5/ 3. Observationally, 



7^ = 1.20±0.06 for clusters with a cool core ( ]De Grandi fc Molendill2002l ). Thus, the heating 
function H st (oc p 3 ~ 7 sA) is similar to the cooling term (oc p 2 A) in equation (J3]). Because of 
this, the balance between the heating and the cooling can be maintained. 



However, the CR streaming is not perfectly locally stable as was noted by lLoewenstein et al 

( I199ll ) (see their § 2). Because of this, the ICM becomes unstable at the end of calcula- 
tions in some models without thermal conduction (e.g. Models LCRe2 and LCRn2). The 
l ocal instability may be related to emission-line filaments, which may be heated by CRs 
(ILoewenstein et al.lll99ll ). Even so, the ability to keep the ICM stable for a long time makes 
CR streaming attractive as a heat source of cluster cores. Moreover, the ability makes it 
easier for the model to attain more stability when it is combined with minor thermal con- 



duction (Models LCRc). Weak tur bulence may also stab i lize the ICM, becau se it conveys 
energy as thermal conduction does (IKim &: Narayanl 120031 ; iFujita et al.l l2004a| ) . 



5. Conclusions 

We have studied heating of cool cores of galaxy clusters by CR streaming. As the source 
of CRs, we considered the central AGN in a cluster. The CRs amplify Alfven waves, with 
which CRs move outward in the cluster. The ICM is heated through dissipation of the waves. 

Using numerical simulations, we found that CR streaming can heat the core for a long 
time after radiative cooling becomes effective without assistance of thermal conduction. 
Development of a strong cooling flow is well prevented. CR streaming can effectively heat 
both high and low temperature clusters. This is because CRs can prevail throughout the 
core and their distribution is insensitive to the change of the ICM. Thus, the entire core is 
heated by the CRs. Minor contribution of thermal conduction makes the ICM even more 
stable. 

There are limitations in our simple models. For example, we assumed spherical sym- 
metry of a cluster. In reality, it is likely that CRs are injected anisotropically. Whether this 
anisotropy is erased or not during the propagation of the CRs in the ICM may depend on 
the geometry of the magnetic fields on which waves and CRs propagate. Moreover, the CR 
injection or the acceleration itself is highly uncertain (S c in equation jH]). In the future, 7-ray 
Observations may reveal CR spectra, which may give us information on CR acceleration in 
clusters. 

We thank the anonymous referee for useful comments. This work was supported by 
KAKENHI (Y. F.: 20540269, 23540308). 
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Fig. 4. — Profiles of the ratios P c /Pg (solid) and Ps/Pg (dotted) at t = 9 Gyr for 
Model LCRO. 
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Fig. 5. — Relative importance of CR streaming (H st ) and collisional heating (if co ii) at t — 
9 Gyr for Model LCRO. 
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Fig. 7.— The evolution of M for Models LCRel, LCRe2, LCRc, and LCRO. 
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Fig. 11. — (a) Temperature and (b) density profiles for Model SCRO. 
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Fig. 13. — CR pressure profiles for Model LCRO. 
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Table 1. Model Parameters 



Model 


Potential 


B (fiG) 


fc 


e 


V 


LCFO 


P 













LCRO 


P 


10 





2.5 x HT 4 


3.1 


LCRel 


P 


10 





1 x 10~ 4 


3.1 


LCRe2 


P 


10 





5 x 10~ 4 


3.1 


LCRc 


P 


10 


0.1 


5 x 10~ 4 


3.1 


LCRnl 


P 


10 





2.5 x 10~ 4 


2.5 


LCRn2 


P 


10 





2.5 x 10~ 4 


3.5 


LCRb 


P 


5 





2.5 x 10~ 4 


3.1 


SCFO 


V 













SCRO 


V 


10 





1 x 10~ 4 


3.1 



